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Abstract. We have performed a first fully 3-D GRMHD simulation with Schwarzschild 
black hole with a free falling corona. The initial simulation results show that a jet is 
created as in the previous simulations using the axisymmetric geometry with the mirror 
symmetry at the equator. However, the time to generate the jet is longer than in the 2-D 
axisymmetric simulations. We expect that due to the additional freedom in the azimuthal 
dimension without axisymmetry with respect to the z axis and reflection symmetry with 
respect to the equatorial plane, the dynamics of jet formation can be modified. Further 
simulations are required for the study of instabilities along the azimuthal direction such 
as accretion-eject instability 



1. Introduction 

Relativistic jets have been observed in active galactic nuclei (AGNs) and micro- 
quasars in our Galaxy, and it is believed t hat they originate from t he regions 
near accreting black holes or neutron stars 
tigate the dynamics of accretion disks and tl 



[ Meier, Koide, fc Uchida] ]. To inves- 
;he associated jet formation, we use 



a full 3-D GRMHD code. One of the most promising models for jet formation is 
the magnetic-acceleration model piandford fc Payne ] . The magnetic-acceleration 



mechanis m has been proposed no t only tor AGIN jets, but also for protostellar 



jets (see [Meier, Koide, h Uchida|]), 



Recently, Koide, Shibata, & Kudoh (1999) jKoide, Shibata, Sz Kudoh ] have 



investigated the dynamics of an accretion disk initially threaded by a uniform 
poloidal magnetic field in a non-rotating corona (either in a state of steady fall 
or in hydrostatic equilibrium) around a non-rotating black hole. The numerical 
results show that the disk loses angular momentum by magnetic braking, then 
falls into the black hole. The infalling motion of the disk, which is faster than 



2 



Nishikawa et al. 



in the non-relativistic case because of the general-relativistic effect below 3rs (rs 
is the Schwarzschild radius), is strongly decelerated at the shock formed by the 
centrifugal force around r = 2rs by the rotation of the disk. Plasmas near the 
shock are accelerated by the JxB force, which forms bipolar relativistic jets. 
Inside this magnetically driven jet, the gradient of gas pressure also generates a jet 
over the shock region (gas-pressure-driven jet). This two-layered jet structure is 
formed both i n a hydrostatic cor ona and in a steady-state falling corona. Koide 
et al. (2000) [[Koide et al 2000|1 have also developed a new GRMHD code in 
Kerr geometry and have found that, with a rapidly rotating (a = 0.95) black- 
hole magnetosphere, the maximum velocity of the jet is 0.9 c and its terminal 
velocity 0.85 c. All of the previous 2-D GRMHD simulations described here were 
made assuming axisymmetry with respect to the z-axis and mirror symmetry 
with respect to the plane z = 0; the axisymmetric assumption suppressed the 
azimuthal instabilities. 



2. 3-D GRMHD Simulations: Equations and Numerical Techniques 



Our basic equations are those of Maxwell for the fields and a set of general- 
relativistic equations representing the plasma, namely the equations of conser- 
vation of mass, momentum, and energy for a single-component conducting fluid 
Weinberg 1972] , Thornc, Price, and Macdonald 1986 , Koide, Shibata, fc Kudoh ], 

Our simulations employ the equations of general-relativistic MHD together 
with the Simplified Total- Variation-Dimi nishing (ST VD) method. The STVD 
method was developed by Davis (1984) [Davis 1984j for studying violen t phe- 
nomena such as shocks (lNishikawa et al. 1997| , [Koide, Shibata, fc Kudoh |. This 



method is as Lax-Wendoroff's method with additional diffusion term. We have 
checked that th is method respects the ene rgy-conservation law and its propaga- 
tion properties [Koide, Shibata, fe Kudoh], 



3. Initial 3-D GRMHD simulations with a Schwarzschild black hole 



In order to investigate how accretion disks near black holes evolve under the in- 



fluence of accret ion instabilities such as the magnetorotational instability [Balbu! 



& Hawley 1995] and accretion-ejection instability (AEI) [[Tagger Pellat, 1999 [ 
the use of a fully 3-D GRMHD is essential. 



3.1. Initial and boundary conditions 

In the assumed initial state, the simulation region is divided into two parts: a 
background corona around a black hole, and an accretion disk (Fig. la). The 
coronal plasma is set in a state of transonic free-fall flow, as in the case of the 
transonic flows with T = 5/3 and H = 1.3; here the sonic point is located at 
r = 1.6rs- The Keplerian disk in the corona is set in the following way. The disk 
region is located at r > ro = 3rs, \cos9\ < 5 = 1/8. Here the density is 100 times 
that of the background corona (Fig. la), while the orbital velocity is relativistic 

and purely azimuthal: v $ = vk = c/[2(r/rs — l)] 1 ' 2 - (Note that this equation 
reduces to the Newtonian Keplerian velocity =y/ GM/r in the non-relativistic 
limit rs/r <C 1). The pressure of both the corona and the disk are assumed equal 
to that of the transonic solution. The initial conditions for the entire plasma 
around the black hole are: p = p$ c + pdis 
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100p ffc 




r > td 
(r < rj) 



and 
or 



\cote\ < 5) 

I cot0 I > 5) 



( , f (0,0, vk) 

(«r,« fl> ^) = | (1^0,0) 



(r > r D 
(r < r D 



and 
or 



| cote | < 5) 
Icot0| > 5) 



(2) 



where we set <5 = 0.125; the smoothing length is 0.3rg. 

In addition, there is a magnetic fi eld crossing the accretion disk perpendicu- 



larly. We set it to the Wald solution fWald 1974 1, which represents the uniform 
magnetic field around a Kerr black hole: B r = BqcosO, Bq = —aBosm.9 (where a 

is the lapse function, a = (1— rg/r) 1 / 2 ). At the inner edge of the accretion disk, the 
proper Alfven velocity is v\ = 0.015c in a typical case with Bq = 0.3-\/ Pqc 2 , where 
the Alfven velocity in the fiducial observer, va = B j\J p + [Fp/ (T — 1) + B 2 ]/c 2 . 
The plasma beta of the corona at r = 3r$ is j3 = p/B 2 = 1.40. The simulation 
is performed in the region l.lrg < r < 20rs, 0<#<7r, < (j) < 2tt with 
100 x 120 x 60 meshes. The effective linear mesh widths at r = l.lrs and at 
r = 20rg are 5.38 x 10 _3 rg and 0.97rs, respectively, while the angular spacings 
along the polar and azimuthal directions are 5.2 x 10~ 2 rad. A radiative boundary 
condition is imposed at r = l.lrg and at r = 20rg. The computations were made 
on an ORIGIN 2000 computer with 0.898 GB internal memory, and they used 
about 47 hours of CPU time for 10000 time steps with 100 x 120 x 60 meshes. 




(b) 




Figure 1. For the fully 3-D simulation,these panels present the time evolution of the proper mass 
density with the magnetic field (B X ,B Z ) in a transonic free-fall (steady-state falling) corona with 
an initially uniform magnetic field, at t — O.Ors (a), t — 39.3rs (b),t = 60.0ts (c), and t = 73.9rs 
(d). The jet is formed around r = 4.5rs as in the 2-D simulation. 



3.2. Simulation results 

Figure 1 shows the evolution of 3-D simulation performed in the region l.lrg < 
r < 20rg, < 9 < it, and < cp < 2tt with 100 x 120 x 60 meshes. The parameters 
used in this simulation are the same as those of the axisymmetric simulations 
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shown in Fig. 6 of | Koide, Shibata, fc Kudoh ]. In this figure, the colored shading 



shows the proper mass density on logarithmic scales; the vector plots show the 
magnetic field. 

The black circle represents the black hole. Figures la presents the initia. 



conditions, which are the same as in the 2-D simulations [ [Koide, Shibata ; ~fc 
Kudohj. At t = 39.2ts (Fig. lb) comparing with Fig. 6c (t = AO. Org) in flKoide 



Shibata, & Kudob] the jet is less generated. At t = 60.0rs as shown in Fig. lc, the 
jet is generated as 2-D simulation at the earlier time (t = 40.0rs). At t = 73.9ts, 
the jet is clearly created around r = 4.5rs, which is shown by the enhanced 
density (Fig. Id). As in the 2-D simulations the jet is generated in a hollowed 
cylindrical form. 

The delay of jet formation seems to be due to reduction of shock formation 
at r = 2rg caused by the additional freedom in the azimuthal direction. Further 
investigation will be reported elsewhere. 

4. Discussion 

Recently, review articles on magnetohydrodynamic production of relativistic jets 
have reported unsolved questions related to jet formation and it s mechanisms 



which det e rmine the veloci t y of jets and time variations of jet flux | Meier, Koide 



fc Uchidal , |Blandford 200l| , |Meier 2002H 



This simulation result is initial and we will perform more simulations and 
investigate effects of the third dimension. This simulation stu dy will be exte nded 
to understand the different (high/soft and low/hard) states [Fender (2001)] and 



different co mbination of accreting rate m = M / M^dd an d the angular momentum 



j = J/Jmax [Meier 2002]. Further results will be reported elsewhere. 
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